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In the framework of the reaction operator approach we calculate and resum the multiple elastic 
scattering of a fast qq system traversing dense nuclear matter. We derive the coUisional broadening 
of the meson's transverse momentum and the distortion of its intrinsic light cone wave function. The 
medium-induced dissociation probability of heavy mesons is shown to be sensitive to the opacity 
of the quark-gluon plasma and the time dependence of its formation and evolution. We solve the 
system of coupled rate equations that describe the competition between the fragmentation of c- 
and 6-quarks and the QGP-induced dissociation of the D- and _B-mesons to evaluate the quenching 
of heavy hadrons in nucleus-nucleus collisions. In contrast to previous results on heavy quark 
modification, this approach predicts suppression of B-mesons comparable to that of D-mesons at 
transverse momenta as low as pr ~ 10 GeV. It allows for an improved description of the large 
attenuation of non-photonic electrons in central Au+Au reactions at RHIC. 
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I. INTRODUCTION 



One of the important experimental signatures of the 
Quark Gluon Plasma (QGP) creation in heavy ion col- 
lisions is the detailed suppression pattern of high trans- 
verse momentum hadrons [l], Q • Jet quenching for light 
mesons, such as 7T, K and 77, at the Relativistic Heavy Ion 
Collider (RHIC) is well explained by radiative energy loss 
calculations [1, 0|- I n contrast, recent data [E IE 0,111 on 
the suppression of single non-photonic electrons cannot 
be reproduced by medium-induced gluon bremsstrahlung 
phenomenology in conjunction with a physically rea- 
sonable set of QGP temperatures and densities d, [Tol ] . 
CoUisional energy loss [ill, EH EH has been proposed 
as a mechanism to reduce the discrepancy between the 
experimental measurements and the perturbative QCD 
(PQCD) estimates of heavy flavor quenching. Still, up 
to transverse momenta pr ~ 10 GeV partonic energy 
loss calculations under-predict the suppression of the 
non-photonic electrons by close to a factor of two [3]. 
The cause of this disagreement is readily identified with 
the small quenching of B-mesons, which dominate the 
high-pT single electron yields. Similar results have been 
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found in a Langevin approach to heavy quark diffu- 
sion where strong heavy quark interactions with D- and 
B-resonances near the phase transition have been as- 
sumed HE EH- 

The model for predicting the quenched hadron spectra 
in A + A collisions, where final-state in-medium effects 
are important, has been very simple: it is assumed that 
the hard jet hadronizes in vacuum, having fully traversed 
the region of dense nuclear matter, L^ GP < 6 fm, and 
lost energy via radiative and coUisional processes. Let 
us examine the validity of this assumption for different 
species of final-state partons and decay hadrons. If a 
parton of non-zero mass ttiq fragments into a massive 
hadron rrih and a secondary light parton, 
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In Eq. (Jlj), P + is the large light cone momentum of the 
parton, z = k+/p+ and |k| ~ A QC d ~ 200 MeV is the 
deviation from collinearity. We can estimate from the 
uncertainty principle for the variable conjugate to the 
non-conserved light cone momentum component Ap~ = 
(p-) f ~ (p-)i-. 
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With Ay + = T form + z form , z form = f3 Q T ioim being the 
formation distance, the formation time reads: 
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Clearly, the higher the energy, or px, and the lighter 
the hadron - the better the perturbative treatment of 
partonic energy loss prior to hadronization will be [l7l |. 
For a px = 10 GeV pion at mid-rapidity Tf orm « 25 fm 
consistent with the jet quenching assumptions 
and the robust description of the nuclear modification 
factor Raa(pt) versus ^/snn and centrality 0, 0, II]- 
In contrast, B- and D-mesons of the same px have for- 
mation times Tf orm ss 0.4 fm and 1.6 fm, respectively, 
<C iy GP . For the purpose of these estimates we used 
z(tt) = 0.7, z(D) = 0.82 and z(B) = 0.92 and the ex- 
act hadron masses in Eq. Therefore, at the finite pr 
range accessible at RHIC and LHC a conceptually dif- 
ferent approach to the description of D- and B-meson 
quenching in A+A collisions is required, when compared 
to light hadrons. In the future in the limit of high trans- 
verse momenta, pt(D) > 10 GeV and pr(B) > 30 GeV, 
improved PQCD description of heavy flavor suppression 
should seek to also incorporate the traditional partonic 
energy loss mechanisms. 

This Letter summarizes key results from our calcula- 
tion of the collisional dissociation of heavy mesons in the 
quark-gluon plasma. In Section II we determine the base- 
line cross sections for heavy flavor production on the ex- 
ample of the identified D° , D + and B + spectra from 
the Tevatron. The light cone wave functions for D- and 
-B-mesons are constructed in Section III. In Section IV 
we use the Gyulassy-Levai-Vitev (GLV) reaction opera- 
tor approach to derive the collisional broadening of the 
propagating qq system and evaluate the survival and dis- 
sociation probabilities of heavy mesons as a function of 
the cumulative transverse momentum transfer from the 
medium. In Section V we solve the system of coupled 
rate equations for the fragmentation into and dissocia- 
tion of heavy mesons. Phenomenological comparison to 
RHIC data on the quenching of non-photonic electrons 
is presented. Our conclusions are given in Section VI. 



II. D- AND B-MESON PRODUCTION IN P+P 
COLLISIONS 

Single inclusive hadron production can be calculated 
to lowest order in the PQCD collinear factorization ap- 
proach as follows: 
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In Eq. Q S = {pa + Pb) 2 is the squared center of 
mass energy of the hadronic collision, x a = p^/pa + , 



1 — 1 — I — 1 — I — 

O B* at s" z = 1 .8 TeV 
□ B* at s" z = 1 .96 TeV 

1=1 PQCD, K = 1 - 1.5 




15 20 
P T [GeV] 



10 15 20 
P T [GeV] 



10 
10 
10 
10 

10' 
10 
10' 
10' 
10 
10 1 



— i — 1 — r 

POCD c-quarks 
PQCD b-quarks 



~\ — 1 — r 

PQCD, r B . 0.02, r D . 0.06 _. 
POCD, r. - 0.07, 1, 




P T [GeV] 



FIG. 1: Perturbative QCD calculation [2l( of the D- and B- 
meson production cross sections at ^J~s = 1.96 TeV, \y\ < 1. 
The top panels show comparison to D° and D + experimen- 
tal measurements at the Tevatron [22l ] and B + data [231. 124^. 
In the bottom panel, the PQCD calculation confronts the 
V^i = 200 GeV non-photonic electron data by PHEHIX at 
RHIC [2^]. The corresponding c- and b-quark cross sections 
are also shown. 



Xb = Pb /pb~ are the light cone momentum fractions 

of the incoming partons and <f>cJN{x a ), <t>b/N(xb) are the 

parton distribution functions [181 ]. In our calculation 

the factorization and renormalization scales are set to 
,2 _ m 2 1 ™2 



/ij = pf, = p\ c + ttlq . The running of the strong cou- 
pling constant is taken to lowest order with V/ = 4 active 
quark flavors and we use m c — 1.3 GeV, nit, = 4.5 GeV. 
Fragmentation functions Dh/ c (z) of heavy quarks into 
mesons, derived in [l||[2(|, are used. At lowest order we 
include Q + g — > Q+g, Q + q—> Q + q and Q + q — > Q + q 
processes, using standard heavy quark PDFs, in addi- 
tion to gluon fusion and light quark annihilation. We 
find that these processes give a dominant contribution to 
the single inclusive D- and B- meson cross section. They 
ensure faster convergence of the perturbation series and 
small next-to-leading order iVfactors. Further details of 
the perturbative calculation are given in (2lj . 

Figure [1] shows comparison of the PQCD single in- 
clusive spectra, Eq. (|4]), to Tevatron measurements at 
= 1.8 GeV and 1.96 GeV [H M, HI. We observe 
that a similar description of the B + cross section can 
be achieved when compared to D° and D + cross sec- 
tions and note that 6-quark fragmentation functions are 
harder (r = 0.07) than the c-quark fragmentation func- 
tions (r = 0.20). Here, r is the parameter that deter- 
mines the hardness of these heavy quark fragmentation 
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functions [l^, H(| and tb/td ~ m c /mf,. We also an- 
ticipate that at px < rrib large corrections will arise to 
the LO perturbative calculation of £?-mesons. However, 
in this region da B j dydpT < da D / dydpT has a smaller 
contribution to the heavy meson production rate. We 
use the same choice of scales and K = 1.5 as the base- 
line for studying nuclear effects on heavy flavor produc- 
tion at the smaller ^snn = 200 GeV at RHIC and the 
larger ^/snn = 5.5 TeV at LHC. We note that if- factors 
cancel in the nuclear modification ratio Raa {pt ) ■ For 
baryons, which constitute ~ 10% of the total c- and b- 
quark decays, we use softer fragmentation functions with 
7"baryon = 1.5r meS cm< In the bottom panel of Fig. [I] results 
for the electrons, 0.5(e + + e~), from the semi-leptonic de- 
cays of heavy flavor, evaluated in LO PQCD in p+p col- 
lisions at RHIC, are shown. Data is from PHENIX 25]. 
Harder fragmentation functions, rg — 0.02, td = 0.06, 
have also been used for direct comparison to alternative 
calculations. 



III. LIGHT CONE WAVE FUNCTIONS 

In order to calculate the effect of the QGP medium 
on a meson traveling through, we need a model for the 
meson state in terms of its partonic degrees of freedom 
(afi,kj). Here x. t = kf / 'P + , J27=i x i = 1 are tne light 
cone momentum fractions and lq, Y17=i^i = are 
the internal parton transverse momenta. To determine 
the typical momenta (kf), we examine the numerical re- 
sults of potential model calculations of D- and i?-meson 
mass spectra and decay widths [26j]. Solving the Dirac 
equation for the relativistic light quark in the potential 
of the heavy quark V = — £/r + br, one can achieve 
a very good description of the lowest-lying and excited 
heavy meson states. The radial wave function of the light 
quark p(r) = r \ip(r)r_ is found to have its maximum at 
clq = 2 — 3 GeV -1 [27]. Fourier transforming a Gaus- 
sian distribution, which features the same maximum of 
the radial wave function p(r), we determine the momen- 
tum distribution: ip(r) ~ e r /( 2a S) — > ip(k) ~ e k a <^ 2 . 
Evaluation of the mean transverse momentum squared is 
straightforward: 

(k 2 ) = N%2nj dcos(d)cos 2 (d) 

x / k 2 kdk e k2a ° = — y , (5) 
Jo 2a 2 

where N 2 ensures the proper normalization of the distri- 
bution in momentum space. Taking do = 2.5 GeV -1 for 
the light parton we obtain (k 2 ) = 0.08 GeV 2 . 

To determine the longitudinal momentum fractions of 
the quarks in a boosted heavy meson we express the 
state in terms of its multi-parton Fock components as 
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FIG. 2: Surface plots of the momentum density \ip(x, Ak/2)\' 2 
for the lowest lying Qq (or Qq) Fock component of the D- and 
B-meson states. 



follows [H S3]: 

x S 2 ^X^kiJ V>(k , (6) 

where the light cone wave functions ^(kj, Xi) of the par- 
tons in the meson are universal. Note that these do not 
depend on the external meson momenta, P + , P. We have 
chosen the following normalization for the single parton 
states: 

(x'P+,k';l|l;k,a;P + ) = (2n) 3 2x S(x - x')6 2 (k - k') . (7) 

In our calculations we take the lowest lying n = 2 Fock 
component. A heavy meson moving in the positive light 
cone direction will then be described by 

|?/>(K, Ak; x, mi, m 2 )| 2 = Norm 2 e i^T^W 2 

x6 2 (K) , (8) 

where the overall large light cone momentum has been in- 
tegrated out. In Eq. |(5J) K = ki + k 2 and Ak = ki — k 2 . 
We assume here that x is the momentum fraction car- 
ried by the heavy quark Q and use m u = m<j = 0.005 
GeV, m c = 1.3 GeV and m b = 4.5 GeV. The light cone 
wave function satisfies the requirement that the max- 
imum in the longitudinal momentum density distribu- 
tion is achieved when the constituent partons of the me- 
son are at the same rapidity or, equivalently, m^/xi — 
m T 2 / x 2 [28;, 29] . Fixing the light quark momentum frac- 
tion squared ((Ak/2) 2 ) to the value obtained in Eq. |(SJ), 
we can determine A = 0.735 GeV and A = 1.055 GeV for 
-D-mesons and i?-mesons, respectively. 

Figure [2] show the momentum density distribution of 
the heavy quark inside the heavy mesons from Eq. ([5]) 
plotted versus the light cone momentum fraction x and 
its half relative transverse momentum Ak/2 = |Ak/2|. 
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The following normalization, J d 2 A\idx \ip(x, Ak)| 2 
I . was used. We note the momentum distribution is 
fairly narrow in Afc/2 and momentum transfers from the 
medium with /j, > 1 GeV 0, i, Q at the initial stages of 
the evolution of the QGP density may easily dissociate 
the heavy mesons if they tend to form early <C L® GP , 
see Eqs. (pj and ([3]). Integrating over Ak we obtain the 
parton distribution function of the heavy quark inside 
the meson: 



4>q/m{x) 



J rf 2 Akd 2 K|^(K,Ak;x,TOi,m 2 )| 2 



Norm 2 47rx(l - x)K 2 e 



-if (l-x) + m^a: 



(9) 



We note that the heavy c- and 6-quark distributions, 
Eq. ([9]), closely resemble in shape the fragmentation func- 
tions D h /q(z) [19|, l20( and peak toward larger values of 
x with increasing heavy quark mass. 



IV. COLLISIONAL DISSOCIATION OF HEAVY 
MESONS 

The GLV reaction operator formalism [3(| HH was 
developed for calculating the induced radiative energy 
loss of hard quarks or gluons when they pass through 
a dense medium. In this approach, the multi-parton 
dynamics is described by a series expansion in \ = 

L QGP OCP 

Jo T &ei(z)p(z)dz = Vf /(A), the mean number of 
interactions that a fast projectile undergoes along its tra- 
jectory. Each interaction is represented by a reaction op- 
erator that summarizes the unitarized basic scattering 
between the propagating system and the medium. The 
summation to all orders in opacity, \i is achieved by a 
recursion of the reaction operator. For the case of col- 
lisional interactions of individual fast partons in dense 
QCD matter, their diffusion in transverse momentum 
space has been derived to leading power [32] and leading 
power corrections (33|. 

The transverse momentum transfer q„ at position n to 
a fast parton when it scatters on the soft constituents of 
the medium is distributed according to the normalized 
differential cross section 
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FIG. 3: Forward cut diagrams representing one collisional in- 
teraction of a qq system with momenta k\ and h,2 propagating 
through dense nuclear medium. 



For further details on the derivation of this formalism, 

see dsn Hi Si. 

The diagrams that are relevant to a single in-medium 
interaction of the quark-antiquark system are shown in 
Fig. [31 We work in terms of the momenta K and Ak, 
defined in the previous Section. If the momentum distri- 
bution of the qq system that has undergone n scatterings 
is oc \M*(K, Ak)M„(K, Ak)|, it can be related to the 
momentum density prior to the last collision as follows: 
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Here e _q "' VK Ak are the momentum shift operators [3l| . 
In Eq. (fT2^) we assume that the overlap between the am- 
plitude and its conjugate varies slowly with q„ and sym- 
metrize the momentum shifts in the last two virtual terms 
around q„ = 0. This allows us to write at the level of 
the momentum distributions (squared amplitudes): 



In Eq. (JTDJ) , jit(n.) = gT(n) is the thermally generated 
Debye screening mass, <J e i{n) rj C^ita 2 ./ jjr(n) with 
Cij = 9/4, 1, 4/9 for gg, qg, qq, respectively, and the mean 
free path A(n) = l/(a e i(n)p(n)). Two momentum trans- 
fers are necessary to build one power of the elastic scatter- 
ing cross section in Eq. (fTOf , allowing for three t — oo on- 
shell cuts in the forward scattering Feynman diagrams, 
see Fig. [3j Momentum flow in the legs of the direct- and 
virtual- interaction terms is constrained as follows: 



Dir.~,5 2 (q„-q: i ) , Vir. 
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|M„(K,Ak)| 2 cx 
2 j g-qn-VAk _ l 



2 f e - q " VK - ij cosh ^ e - q " v ^ 
M„_!(K,Ak)| 2 . (13) 



In Eq. (|13[) we also used the symmetry of the momentum 
transfer distribution, Eq. (fTO]) . relative to the transfor- 
mation q„ — ► — q„. The basic step in Eq. (|T5j) allows 
us to resum the interactions to all orders in opacity and 
relate the final momentum distribution of the evolved qq 
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system to the one in Eq. (JSj> : 

i^.(K,Ak)p = E^/n^ 



1 da, 



e q " VK - 1 ) cosh f q„- V Ak 
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An approximate closed form for Eq. (fl4|) can be ob- 
tained by Fourier transforming to the impact parameter 
space (B,b) conjugate to (K, Ak): 



|^/(B,b)| 2 = |MB,b)| 2 exp 2 X / d 2 q 
x ( (e iq B - 1) cos (q • b) + (e lq b - l) ) 



1 da. 



a e i d 2 q 
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A closed form for the integral in Eq. ([15]) does not exist 
even for simple forms of the differential elastic scattering 
cross section. The cosine term couples the broadening 
of the total momentum K of the quark- antiquark system 
to the distortion of the light cone wave function in Ak. 
We can still calculate the final momentum distribution 
by considering only the leading effect of the coupling: 
cos(q • b) = 1 — (q ■ b) 2 /2 + ■ • • . Next, we perform the 
integrals over the Yukawa potential. With q = |q|, 
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The b — > B contribution without interference yields a 
result similar to Eq. (fT!)]) and the second term in the 
expansion of the cosine gives 
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In Eq. (fl~7|) we have taken B and b in the same direction 
and chosen 1 / B as the upper limit of the q integral to 
obtain an estimate for the upper limit of the interference 
contribution. The key to evaluating the average K- and 
Ak-broadening of the partons is the small bfi expansion 
in Eqs. (fTC)) and (TIT)) : 
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Note that the leading correction that arises from Eq. (TTT)) . 
3& 2 /^ 2 /16, is small. Keeping terms cx XM 2 £> where £ = 
ln(2e~~* E /(b^)) + 1/2 - 3/16 > 0(1), and incorporating 
only the leading effect of the coupling term in Eq. (TTT|) 
we find: 
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We treat £ as approximately constant when compared to 
the power behavior of b 2 and B 2 when Fourier transform- 
ing Eq. (|19[) back to momentum space: 
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In providing a physics interpretation to Eq. (f^TTj) , we note 
that the first part of our result represent the broaden- 
ing in the momentum distribution of the meson itself, 
(K 2 ) = 4xm 2 £- It is twice the size of the broadening for 
an individual parton [33| . The second part, critically im- 
portant for this work, represents the distortion in the in- 
trinsic momentum distribution of the quarks, which leads 
to the meson decay. 

We first integrate out the distribution in K. The 
small acoplanarity that may arise in heavy meson or 
non-photonic electron triggered correlations [U, [34| is 
neglected in this work. Since our wave functions are real, 
they are given by the square root of the Ak and x momen- 
tum distributions, Eqs. ([5]) and PU|) . In our calculation 
the initial-state ip (Ak,x) = -0A/(Ak, x) represents the 
inclusive D- or S-mesons, respectively. The final-state 
ipf{Ak,x) = a^i/(Ak, x) + (1 — a)^(Ak,i) denotes 
a superposition of the meson and a dissociated qq pair. 
The survival probability is given by a 2 if the light cone 
momentum distributions are normalized to unity. We 
readily obtain: 
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Eq. (PT|) is one of the main theoretical results derived in 
this Letter. Although the integral over the heavy quark 
light cone momentum fraction x has to be taken nu- 
merically, direct comparison of the integrand to Eq. (J5|) 
shows that the meson survival probability P a (XM s) — !• 



Equality is reached when x/i 2 £ = 0, i.e. in the absence 
of interactions in the medium. We finally note that typ- 
ical values of £ can be estimated from the requirement 
bfi -c 1: for bp = 0.25 — 0.1 wc find £ ~ 2 — 3, respec- 
tively. 



V. HEAVY MESON SUPPRESSION 
PHENOMENOLOGY 

We are now ready to calculate the suppression of heavy 
hadrons, H(c) and H(b), from collisional interactions in 
the QGP. The fragmentation time, averaged over the 
final-state mesons and baryons is calculated as follows: 
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In Eq. (|2"2")) the fragmentation functions are normalized to 
the fragmentation fractions for Q — > H(Q), which in the 
QCD factorization approach are assumed to be universal, 

Jp 1 dzD H,/Q(z) = fffiiQ), E l /ffi(Q) = !• The dissocia- 
tion time, on the other hand, is externally driven by the 
dynamics of the evolving bulk medium and calculated 
by taking the logarithmic derivative of the dissociation 
probability, 



Pd(PT,m Q ,t) = 1 - P s (p T ,m Q ,t) , 



derived in Eq. (|2"Tj) . as follows: 
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We use the same initial soft gluon rapidity density 
dN 9 /dy as in the calculation of the 7r° quenching [4| in 
central Au+Au and Cu+Cu collisions at RHIC and cen- 
tral Pb+Pb collisions at LHC and the cumulative mo- 
mentum transfer is given by 



Aq to 



(25) 



In Eq. (|2"5"]) to < t and to = 0.6 fm, consistent with 
QGP formation times used in hydrodynamic simulations 
of bulk observables [35], and (3q = dz/dt. Let us denote 
by 



f Q (PT,t) = 
f H (pr,t) = 



daQ{t) 

dyd 2 pT 
da H (t) 
dyd 2 p T 



f Q (p T ,t = 0) = 
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dyd 2 px 
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(26) 
(27) 



the double differential cross sections for the heavy quarks 
and hadrons (mesons+baryons) . Initial conditions are 
also specified above, in particular the heavy quark dis- 
tribution is given by the perturbative c- and 6-quark jet 



0.75 - 



Q. 

~~5 0.5 

< 



0.25 




0.75 - 



Q. 

~~$ 0.5 

< 

cc 



0.25 - 



. 200 GeV 

B+D mesons+baryons 
B mesons+baryons 
D mesons+baryons 



Central Au+Au at RHIC (dN°/dy = 1 1 75) 
J I I 



Central Pb+Pb at LHC (dN°/dy . 2000) 
_1 I I I I 




Central Cu+Cu at RHIC (dN°/dy . 350) 

, L 




Central Pb+Pb at LHC (dN J /dy . 3500) 
I I I I I 



10 20 
P T [GeV] 



10 20 30 

P T [GeV] 



FIG. 4: Suppression of D- and jB-meson production via col- 
lisional dissociation in the QGP. Top panels show numerical 
results for Raa(pt) in central Au+Au and Cu+Cu collisions 
at RHIC for gluon rapidity densities dN 9 /dy = 1175 and 350, 
respectively Q. Bottom panel shows predictions for central 
Pb+Pb collisions at LHC for dN 3 /dy = 2100 and 3500. 



cross section. The fragmentation fraction of &-quarks into 
i? c -mesons is very small. In our work it is neglected and 
the rate equations that describe the competition between 
b- and c-quark fragmentation and D- and jB-meson dis- 
sociation decouple for different heavy quark flavors. In- 
cluding the loss and gain terms we obtain: 



dtf Q (PT,t) = - 



1 



-f Q ( P T,t) 



1 



(Tdiss(pT/x,t)) Jo X 



T£aim(PT,t)) ' 
1 I 

dx—(f> Q/H (x)f H (pT/x,t) , (28) 



dtf H {pT,t) 



1 



(T diss (p T ,t)) 



f H (PT,t) 



(TiatmipT I Z,t)) Jo z 



dz-D H/Q (z)fQ( PT /z,t) .(29) 



In Eqs. (|28|) and (I29| z and x are typical fragmenta- 
tion and dissociation momentum fractions and we have 
checked that in the absence of a medium, Tdi ss (pT, t) — > 
co, we recover the PQCD spectrum of heavy hadrons 
from vacuum jet fragmentation. 

We solve the rate equations, Eqs. (|2"5]) and (I29I1 . nu- 
merically using the initial conditions, Eqs. (f2"6"| and (|27[) . 
and obtain the final spectra of heavy hadrons at t > 
Ly GP ,Tf orm and Tdi ss - The nuclear modification factor, 
Raa{pt)i which arises from the collisional dissociation 
of £>-and £?-mesons, is shown in Fig. 0] for physical sit- 
uations expected to be prevalent at RHIC and LHC. In 
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spite of the inherently different physics mechanisms that 
drive the suppression of light and heavy hadrons, com- 
parison to the quenching of the 7r° [3] can be done if the 
same model for the QGP properties, such as the gluon 
rapidity density dN 9 /dy, is used. The theoretical uncer- 
tainty comes from varying the parameter £ in its natural 
range and /io and Ao are set by the overall gluon multi- 
plicity and the QGP formation time, To, using thermody- 
namic arguments. The top panels show results for central 
Au+Au and Cu+Cu collisions at the maximum RHIC 
energy ^snn — 200 GeV. The bottom panels present 
Raa(pt) for two different dN 9 /dy extrapolations in cen- 
tral Pb+Pb collisions at LHC. 

One important feature of this approach is the sensi- 
tivity to the build-up of the QGP density. If To << 
Tform and p(t) decreases rapidly, hadrons would not have 
formed at times when the dissociation mechanism is most 
efficient. This leads to reduced suppression, which can be 
seen in the high-p^ behavior of the D-meson Raa and 
understood when one recalls that Tf orm oc px- Since the 
distortion of the light cone wave functions reflects the 
accumulated squared transverse momentum transfer, the 
reduced suppression at px < rriQ arises from the velocity 
factor (3q in Eq. (|25|) . We observe in Fig. [4] that with 
initial formation time To = 0.6 fin, consistent with the 
one used in hydrodynamic simulations to describe bulk 
QGP observables [35|], the moderate- and high-p^ sup- 
pression of heavy hadron production is comparable to 
that of lighter hadrons. We have checked that for QGP 
formation time To = 0.2 fm the px > 5 GeV heavy quark 
quenching is reduced by a factor of ~ 1.5. Finally, one 
should note that dissociation/fragmentation both emu- 
late energy loss by shifting the quarks/hadrons to lower 
transverse momenta. For example, B- and D-meson at- 
tenuation is sensitive to the partonic slope. For this rea- 
son the suppression at LHC is found to be comparable or 
slightly smaller than the one calculated at RHIC in spite 
of the larger QGP densities and temperatures. 

Contrary to calculations that emphasize radiative and 
collisional heavy quark energy loss [t| EH EH EH , QGP- 
induced dissociation predicts S-meson suppression com- 
parable to or larger than that of D-mesons at transverse 
momenta as low as px ~ 10 GeV, see Fig. [4] This is due 
to the significantly smaller formation times for H (b) rel- 
ative to H (c). Thus, each fragmentation/dissociation cy- 
cle proceeds at a much faster rate for fr-quarks/B-mesons. 
This is an example where the large mass facilitates the 
hadron suppression mechanism. The observable effect of 
this faster rate in the quenching of the final hadron distri- 
butions is amplified by the significance of the early t ~ To 
hot and dense stage in the dynamical evolution of the 
QGP. Our simulations show that the rates in Eqs. (T2"51) 
and (f2T)|) may play a more important role than the ef- 
fective fractional energy losses e = AE/E w 1 — z and 
e w 1 — x, which are larger for c-quarks/D-mesons. 

Inclusive non-photonic electron data for the most 
central Au+Au collisions at ^/Jnn — 200 GeV from 
PHENIX H| and STAR are compared to the current 
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FIG. 5: Suppression of inclusive non-photonic electrons from 
D- and B-meson spectra softened by collisional dissociation 
in central Au+Au collisions. Data on non-photonic electron 
quenching from PHENIX and STAR is also shown. 



calculation in Fig. [5] We used the PYTHIA event gen- 
erator [36J to simulate the (e + + e _ )/2 spectra coming 
from the primary decays of the D- and B-mesons in our 
baseline p+p and QGP-modified A+A results. We find 
that the single electron suppression can be as large as a 
factor of four and approximates well the quenching ex- 
tracted from available data 0, H, 0, [|j|. We emphasize 
that such agreement between theory and experiment is 
not achieved at the cost of neglecting the contribution of 
the £>-mesons to the non-photonic electron spectra. 



VI. CONCLUSIONS 

It has been recently suggested that collisional dissoci- 
ation of heavy quarkonia in the quark-gluon plasma [37j 
may be a possible explanation for the suppression of their 
production rate in nucleus-nucleus collisions. It is thus 
surprising that until now a similar physics mechanism has 
not been considered for open heavy flavor. In this Letter, 
we investigated the perturbative QCD dynamics of open 
charm and beauty production and, in the framework of 
the reaction operator (GLV) approach [3(| HH HH, HH 
extended to composite qq systems, derived the medium- 
induced dissociation probability for heavy D- and B- 
mesons traversing dense nuclear matter. We showed that 
the effective energy loss, which arises from the sequen- 
tial fragmentation and dissociation of heavy quarks and 
mesons, is sensitive to the interplay between the forma- 
tion times of the hadrons and the QGP and the detailed 
expansion dynamics of hot nuclear matter. The pro- 
posed new attenuation mechanism, which stems from the 
short formation times of D- and B-mesons and under- 
lies the suppression of the inclusive non-photonic elec- 
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trons 0, [f| @, [1] in nucleus- nucleus collisions at RHIC, 
was found to be compatible with the measured large, fac- 
tor of four to five, quenching for heavy flavor. 

Previous studies, based on radiative and collisional 
parton energy loss 0, [l(| Hi and heavy quark diffu- 
sion [lH ] under-predict the suppression of non-photonic 
electrons in central Au+Au collisions due to the small 
fo-quark quenching. A natural consequence of the ap- 
proach, presented in this Letter as a viable alternative 
to existing calculations, is that i?-mesons are attenuated 
as much as L>-mesons at transverse momenta as low as 
Pt ~ 10 GeV. While we anticipate that a comparable de- 
scription of the attenuation of the non-photonic electrons 
may be achieved when partonic energy loss is combined 
with quark-resonance interactions near the QCD phase 
transition in Langevin transport simulations, the hierar- 
chy Raa(H(c)) <§; Raa(H '(b)) will not be changed in the 
accessible transverse momentum range [38[ . We conclude 



that robust experimental determination of the dominant 
mechanism for in-medium modification of open heavy fla- 
vor would require direct and separate measurements of 
the B- and D-meson Raa distributions versus px and 
centrality in collisions of heavy nuclei. 
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